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1. INTRODUCTION AND SUMMARY 


This report summarizes work performed on JPL Contract 956999, 
"Development of Single-Crystal Silicon Dendritic Ribbon and High- 
Efficiency Solar Cells.” The contract was initiated September 12, 1984 
and was completed on December 31, 1984. Work reported herein was 
completed under a "Partial Order Transfer" from the Westlnghouse 
Advanced Energy Systems Division to the Westlnghouse Research and 
Development Center. The bulk of this work was performed at the R&D 
Center, and this report represents work performed both at the R&D Center 
and at AESD. 

Work performed under this contract was, in fact, a direct 
extension of work performed under JPL Contract 955843. This effort has 
continued into 1985 under yet another contract (957207). The long-range 
objectives of all these contracts are the same: to conduct development 
efforts to achieve initial deliverables of a JPL program aimed at 
demonstrating that the silicon dendritic web technology is ready for 
commercial use by the end of 1986. A commercial readiness goal involves 
improvements to crystal growth furnace throughput to demonstrate an area 
growth rate of greater than 15 cm /min while simultaneously growing 
10 meters or more of ribbon under conditions of continuous melt 
replenishment. "Continuous" means that the silicon melt is being 
replenished at the same rate that it is being consumed by ribbon growth 
so that the melt level remains constant. 

Specific tasks of this contract No. 956999 Included computer 
thermal modeling required to define high-speed, low-stress, continuous 
growth configurations; a study of convective effects in the molten 
silicon and growth furnace cover gas; furnace component modifications; 
web quality assessments; and experimental growth activities to 
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demonstrate progress. A specific milestone identified in this contract 
involved demonstration of an area growth rate greater than 10 cm^/min 
while simultaneously growing 10 meters or more of ribbon under 
conditions of continuous melt replenishment by the end of 1984, 

Overall progress toward achieving increased crystal lengths in 
this as well as previous related contracts dating from 1980 is presented 
in Figure 1. Crystal lengths of up to 10.9 meters were achieved prior 
to the end of 1984 during growth with no replenishment. Replenished 
growth configurations have produced crystals up to 7 meters — the 1984 
goal being 10 meters. 

Figure 2 further elaborates on progress made in the web growth 
program in 1984 (this contract covering the last three-month portion of 
the year) toward goals in continuously replenished crystal growth. The 
maximum crystal length grown in a melt replenishment mode prior to March 
1984 was 2 meters; the crystal required three hours of growth time. 
Through the course of this program, crystal lengths have more than 
tripled; growth times, before interruption of growth due to the onset of 
thermal instability, have increased by almost a factor of four. 

The 1984 web growth goals of crystal length and time are also 
presented in Figure 2. A comparison will show that stability 
improvements to increased uninterrupted growth times have exceeded the 
goals. However, crystal pull speeds have not increased as was hoped, 
resulting in shorter crystals than projected. The pull speed is limited 
by stresses in the growing crystal that can themselves interrupt 
growth. As a result of work completed on this contract, the primary 
focus of the 1985 Web Growth Development Program is being directed 
toward the solution of the stress problem. 

The overall status of the Web Growth Development Program at the 
end of 1984 is presented in Table 1. 
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Figure !• Progress made toward growth of long crystals. 
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Figure 2. Continuously replenished web growth. 
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Table 1. Web Growth Development Status — December 1984 


Primary Thermal Variable Identified and Fix Verified 

Significantly Improved Thermal Isolation of Growth and Melt Regions 
Achieved 

2 

6,8 m Long Ribbon Grown at 7 cm /min Rate with Replenishment 

Replenished Web Production Equivalent to Pre-Pilot Line Furnaces 

Computer Based Data System Being Used to Assist in Parameter 
Correlations 

10 ra Long Crystal Grown 

2 

Additional Stress Reduction Work Required to Achieve 10 cm /min Rate 
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2, TECHNICAL PROGRESS 


2.1 Thermal Stress and Buckling Modeling 

Modeling of thermal stresses and buckling of dendritic web 

crystals during growth operations is required to define growth 

configurations capable of continuous production of high-quality ribbon 

(having suitable thickness and low stress) at increased area growth 

rates. In particular, the ultimate goal of this contract was to 

2 

demonstrate an area growth rate of greater than 10 cm /min while 
simultaneously growing 10 meters or more of ribbon under conditions of 
continuous melt replenishment. 

Experimental temperature data from elements of a state-of-the- 
art growth configuration (labeled J460) were used as input for initial 
modeling work. The experiments were performed prior to initiation of 
this contract. Previous modeling of the J460 growth configuration has 
been performed using calculated element temperatures. The principal 
differences were that the newer cases had a much colder second lid (L2) 
and a somewhat hotter shield stack than in previous models. Results of 
the initial calculations performed using experimental temperature data 
indicated that the growth velocity increased by about 16% (compared to 
previous calculations) but that interface stress had concurrently 
increased about 40%. The far stress peak, a principal source of 
buckling, decreased by about 40%. While these results indicated the 
sensitivity of growth parameters to lid temperature changes, it was felt 
to be too drastic a change in L2 temperature. Accordingly, plans were 
made to check the temperature data for the second lid using optical 
pyrometry. 

Based on the effect of the hotter stack temperature, a modifica- 
tion of the J460 configuration was modeled to include a vertical thermal 


6 


element on the top of the normal shield stack. The calculated tempera- 
ture profile for the modified configuration indicated that the far 
stress peak should be almost eliminated. Stress calculations were 
initiated; set-up of an experimental furnace to check the results was 
also initiated. 

Difficulties were encountered in performing a stress or buckling 
configuration for the modified configuration as a result of operating 
system changes in the computers at the Westinghouse Nuclear Center where 
the WECAN code is run. To allay this problem, the temperature calcula- 
tion code was transferred from the Univac 1180 (at the Westinghouse R&D 
Center) to the Cray-1 (at the Westinghouse Nuclear Center), allowing 
direct interface with the CDC7600 that runs the WECAN stress and 
buckling programs. 

In addition, stress modeling calculations were made to investi- 
gate the effect of an asymmetric variation of temperature across the 
width of the ribbon and thus complete the investigation of lateral 
temperature variations in web crystal growth. Two cases were considered. 

In the first case, a linear variation in temperature was imposed 
across the width of the crystal. In the second case, a cubic variation 
in temperature was imposed across the width. In both cases, the 
absolute temperature difference between the centerline and the dendrites 
was ten degrees as in the previous cases. The temperature profile was 
assumed to be: 

T(x,y) = T(x)g(y) (1) 

where T(x) is expressed mathematically by: 

T(x) = (A/B^)exp(-Bx) + Cx + D. (2) 

In these expressions. A, B, C, and D are constants, x is the position 
along the length of the crystal measured from the solid-liquid inter- 
face, and y is measured across the ribbon width from the centerline. 
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The y-variation of temperature functions used were: 


g(y) = 1 + .00593 (y/w) (3) 

and 

g(y) = [1 + .00593(y/2)^]/(1.0015). (4) 

In these expressions, w is the ribbon half-width, 

A significant difference between the calculations performed here 
and any previous ones, at least from the point of computation time, was 
the necessity of using a full-width mesh as shown in Figure 3. The 
extraction of the first and second eigenvalues (buckling modes) 
increased the computer time approximately four-fold over the earlier 
cases. The stresses in the case of the isotherms with a linear 
y-variation were very similar to those calculated previously in models 
in which no y-variations in temperature were considered. This should 
not be surprising since the linear y-variation in temperature would not 
change the Laplacian of the temperature. There was some change in the 
stress pattern for the cubic case with the stress Increasing on one side 
and decreasing on the other. Again, this effect can be understood in 
terms of the effect of the y-variation of the temperature. 

The buckling results are somewhat more surprising. The eigen- 
values for both symmetric and antisymmetric buckling (bending and 
twisting) are given in Table 2. 

Considering that the numerical algorithm currently used for 
extracting the eigenvalues is different from that used when the flat 
isotherm case was analyzed, the small variations in the numerical values 
are insignificant. The lack of any effect from the linear y-variation 
is reasonable since the stress profiles are not affected as noted. The 
insensitivity of buckling to the cubic y-variation is more surprising 
and might be attributed to the lack of any significant net change in the 
stress across the width of the ribbon. 
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Figure 3. Full width WECAN mesh for buckling analysis. 
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Table 2* Buckling Eigenvalues 


Case 


Eigenvalue 


Symmetric 

Antisymmetric 

Flat Isotherms 

1.467 

1.597 

Linear slope in y 

1.448 

1.568 

Cubic in y 

1.445 

1.565 


These results have an important implication for the actual 
growth of the web ribbons in that the asymmetry of the temperature 
profile across the width of the ribbon has negligible effect compared to 
the temperature profile along the length of the ribbon, and for any 
positive or negative net curvature of the isotherms that may exist# 

Analytical work was completed on this contract using a "first- 
order** approach which involved use of separate models to describe 
various aspects of the dendritic web growth system. Aspects considered 
involve temperature distribution in the web crystal during growth, stress 
in the web crystal resulting from these temperature distributions, fluid 
flow in the crucible containing the molten silicon, and the temperature 
distribution in the molybdenum susceptor that heats the silicon. 

Figure 4 is a logic diagram that shows how the web temperature 
and stress analysis are coupled in the **f irst-order modeling** technique. 

Table 3 summarizes the most significant results achieved in 1984 
through the use of the first-order modeling. Conclusions concerning the 
continued use of the first-order approach are summarized in Table 4. 

Based on these results and conclusions, a more comprehensive 
approach has been identified for use in future programs. This approach 
would integrate existing models so that a single model could be used for 
analyzing the web growth system (with the exception of the 
stress/buckling calculations). The advantage would be that the system 
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temperature 
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Figure 4. Logic Diagram of First Order Modeling Approach. 
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Table 3. Results of First Order Modeling 


!• Hardware developed on basis of modeling: 

a. Grows web crystals 5 to 6 cm wide 

b. Grows web crystals at 2 to 3 cm/min 

[But not simultaneously on a long-term basis,] 

2. Model predictions validated when real system temperatures are used 
for boundary conditions, 

3. Model has predicted a design for simultaneous growth at W > 5 cm and 
V > 2 cm/min, 

4, Highest speed growth configurations can produce crystals with large 
residual stress, 

[Residual stress generation is not included in model, but has not 
been a problem to date.] 

5, Simultaneous satisfaction of the required boundary conditions for 
web temperature and melt temperature can be complex. 


Table 4, Calculations Based on First Order Modeling 


1. Interactions of system elements is becoming more important with more 
advanced systems, 

2. A more comprehensive analytical model is required for guidance in 
the integrated design of growth systems. 

3. The significance of residual stress is not yet determined. 


interactions could be readily evaluated and the design and evaluation 
process speeded up. Additionally, the more complete model would allow 
calculation of the effects of system changes that are presently 
accessible only by time-consuming experimental iterations, A logic 
diagram of the advanced modeling concept is shown in Figure 5, 

The three-step procedure used in this program to analyze stress 
and buckling in growing dendritic web is outlined in Table 5. The 
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INTERPRET 

ELASTIC 

MODEL 
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PLASTICITY 
IN STRESS 
MODEL 

T — 


Figure 5 


DESIGN/MODEL 

Idealized 

Configuration 


I MODEL/PREDICT TOTAL 
SYSTEM INTERACTIONS 


BUILD AND TEST 



Similar to First 
Order modeling 

Hay use finite 
element veb teo^> . 
erature model 


3-D System Model 
for temperature 

Fluid flow In melt 


Measure system 
ten^eratures 


Measure residual 
stress In web; does 
It exceed acceptable 
values 


Logic diagram of advanced modeling approach. 
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Table 5. Three-Step Procedure for Analyzing Silicon Ribbon 


• Step 1 — Determining Temperature Profile 

• Step 2 — Performing Thermal Stress Analysis Using Temperature 

Profile Obtained From Step 1 

[K]{A} = {F} 

{a} = [E][D]{A> 

• Step 3 — Performing Buckling Analysis 

[K]{A} = X[Kg]{A} 

X < 1 buckled 
X > 1 not buckled 


stress analysis is completed using the Westinghouse "WECAN" computer 
model. Basic capabilities of the model are summarized in Table 6. 
Figures 6 and 7 show the two- and three-dimensional element mesh 
capabilities, respectively, of the WECAN code. 

Effects of asymmetric temperature variations in the y-direction 
(across the width) of the growing web were analyzed using the WECAN 
model. The two temperature variations imposed for this analysis are 
listed in Table 7. Computer-printed stress contours associated with the 
linear temperature variation are shown in Figure 8. These results are 
essentially identical to stress contours calculated in a previous 
program using flat isotherms. 

Stress contours calculated using the asymmetric cubic 
temperature variation are shown in Figure 9. The asymmetry of the cubic 
contours is evident in the figure. Surprisingly, the critical buckling 
conditions were unaffected by the y-direction temperature profiles 
imposed across the width of the web. 

The calculated ribbon shapes for the symmetric and asymmetric 
temperature variations are shown in Figures 10 and 11, respectively. 
Again, critical buckling shapes are essentially identical to those for 
the flat isotherm case. 
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Table 6. WECAN Analysis Capabilities 


o Static Analysis 
Elastic 

Inelastic (Plasticity and Creep) 

Large Deformation Prebuckling Analysis 

e Thermal and Flow Analysis 

Steady State and Transient Heat Conduction Analysis 
Steady State Flow Analysis 

Other Field Problems (Electromagnetic, etc.) 
o Dynamic Analysis 
o Harmonic Response Analysis 
• Nonlinear Transient Dynamic Analysis 
o Linear Buckling Analysis 



Figure 6. WECAN variable spacing finite element mesh. 


15 





Figure 7. WECAN three-dimensional finite element mesh* 

2.2 Liquid Convection and Heat-Transfer Modeling 

The purpose of this task was to develop a quantitative 
understanding of convection and heat transfer in the liquid pool of 
silicon during dendritic web growth, with the ultimate objective of 
aiding both process optimization and the evolution of an optimal design. 
Work on this task was performed under subcontract to Cambridge Materials 
Modeling Group (CMMG) directed by Prof. Julian Szekely. Most of the 
work was performed by N. El-Kaddah of CMMG. 
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Table 7. Y-Direction Web Temperature Variation 


T(X,Y) - TM (X) F(Y) 

WHERE 

TM(X) = * EXP (-BX) + CX + D 

F(Y) = I + Cl jt (Y/W) + C3 / (Y/W) 

CASE 1 - LINEAR VARIATION (C3=0) 

A = 3000, B = 5, C * -86.5 AND D = 1565 

THIS IMPLIES THAT TM(0) = 1685 AND TM(IO) = 700. 

Cl = 0.00593 WHICH IMPLIES THAT T(0, 1.975) - T(0, 0) = 10 
CASE 2 - CUBIC VARIATION (C1=0) 

A, B, C AND D ARE SAME AS CASE 1. 

C3 = 0.0593 WHICH IMPLIES THAT T(0, 1.975) - T(0, 0) = 10 


The domain of interest of this work is sketched in Figure 12. 

The walls of the crucible are heated while a thin film of silicon 
monocrystal (dendritic web) is being withdrawn, as also Illustrated in 
the sketch. 

The objective of this analysis was to define convection within 
the liquid pool. This convection plays an important role in affecting 
the temperature field and its stability in the vicinity of the dendrites 
at the leading edge of the solid film, as the crystal is being withdrawn. 
This analysis was intended to lead to an understanding of techniques 
required to increase both the film withdrawal speed and the film width 
while minimizing convection in the molten silicon. 

In the initial phase of the convection modeling work it was 
suggested that the most critical temperature distribution was probably 
in the cross section of the melt normal to the plane of the web crystal. 
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Figure 11. Calculated ribbon shape for asymmetric temperature variation 



and a hypothetical, but reasonable, temperature field was selected for 
this region. Using this field as the input data, Szekely and El-Kaddah 
calculated the fluid flow in the silicon melt. A typical example is 
shown in Figure 13, where the temperature field is shown at the top of 
the figure. In the particular case illustrated, a melt depth of 0.9 cm 
was assumed. A similar calculation was made for a melt depth of 0.7 
cm. The contrast between the two cases was striking. The 0.2 cm 
reduction in melt depth resulted in a decrease in the maximum flow 
velocity by a factor of two. This change is in qualitative agreement 
with the experimental observation that reducing the melt height by a 
millimeter or so greatly improves the stability of web growth. 

Parametric variations were then studied to provide a basis for 
possible future crucible designs. Parameters selected for evaluation 
were crucible width, length, and temperature asymmetry. In addition, a 
different temperature profile, calculated from a two-dimensional WECAN 
model of the growth system, was provided as input for evaluation. 
Although this temperature field was qualitatively similar to the one 
used for the initial flow calculations, it was quantitatively different 
such that resulting flow patterns would give some indication of the 
sensitivity of the system to the details of the temperature field. 

The two-dimensional mesh shown in Figure 14 was used to 
calculate the temperature field in the whole furnace system (although 
the melt was, of course, the region of particular interest in this 
analysis). For simplicity, radiation connections between elements are 
not shown in the figure. Input data for the initial exercise of the 
model was the measured temperature profile shown in Figure 15. The 
resulting melt temperature distribution is shown in Figure 16. 

The second case run with the model was a parametric variation in 
which the emissivity of the silicon melt and the molybdenum surface 
elements were increased by 50%, i.e., from 0.2 to 0.3. The same 
temperature profile (Figure 15) was used as the input data for the 
second run. The resulting melt temperature distribution is shown in 
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Figure 13. 


B. Calculated Velocity Field for 
a 9 im Melt Deep Silicon Melt 


Calculated temperature and velocity fields in a silicon 
melt (from Szekely and El-Kaddah). 
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Figure 15. Temperature profile used as input data for finite 
element analysis. 
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Figure 16. Isotherms in silicon melt resulting from temperature input 
of Figure 15 when e = 0.2. 
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Figure 17. Isotherms in silicon melt resulting from temperature input 
of Figure 15 when e = 0.3. 
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Figure 17. The temperature gradients increased approximately in 
proportion to the emissivity change. Based on experimental melt 
temperature data, it appeared that an intermediate emissivity would be 
more appropriate to the physical situation. 

An additional case was run using the two-dimensional model, but 
changing the power input distribution. In the first two cases, the 
power input varied along the wall of the susceptor in such a way as to 
create the specified temperature distribution. In the third case, a 
uniform power input was assumed along the susceptor wall and adjusted in 
magnitude to give the proper temperature at the centerline. Some 
changes occurred in the temperature distribution relative to the first 
case, but they were less dramatic than those resulting from the change 
in emissivity. 

Finally, a parametric study on the effect of crucible length and 
temperature on the convective flow in the silicon melt was performed by 
Profs. Szekely and El-fCaddah. Crucible length is measured in the x- 
direction of the sketch shown previously in Figure 12. The various 
cases examined in this study are shown schematically in Figure 18. The 
first group of cases was formulated to show the effect of the crucible 
dimensions; the second group was formulated to investigate temperature 
effects; the third case Illustrates the effect of front/back temperature 
asymmetry in the melt (across the width of the cruible); and the fourth 
case is a representation of the side-to-side temperature profile in the 
crucible. The centerlines shown in Figure 18 are the positions where 
the seed first contacts the molten silicon (theoretically) during the 
initiation of crystal growth. 

The results of these calculations can be summarized as follows: 

1. For a given temperature difference between the center and 
front (or back) of the crucible, the maximum fluid 
velocity decreases with Increasing crucible length. 

2. For a given crucible length, the maximum fluid velocity 
increases with increasing temperature difference. 
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Figure 18. Dimensions and temperatures used for parametric analysis of 
flow in silicon melt. 
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These two results can be combined to conclude that the fluid flow 
depends on the temperature gradient. Furthermore, the results indicate 
that it is the front-to-back temperature gradient rather than the 
vertical temperature gradient that dominates convective fluid flow in 
the molten silicon. 

The results in case III, the asymmetric gradient, were of 
particular interest. The two-degree asymmetry resulted in a shift in 
the temperature centerline of about 0.7 mm while shifting the center of 
the flow pattern by about 5 mm. As a result, a web crystal growing in 
the center of the melt would be only slightly away from the proper 
thermal position but would be in a location of high fluid flow, where 
deleterious thermal fluctuations could occur. Thus it is important to 
obtain a good front-back thermal balance in web growth systems. 

In case IV, the approximation to the side-to-side temperature 
distribution, the results were not so dramatic but still of considerable 
interest. As would be expected, when the flat temperature profile is 
induced over the central region of the crucible, there is essentially no 
lateral flow in this region. This lack of lateral flow gives more 
veracity to the flow patterns calculated in Cases I through III, since 
those flow patterns would not be distorted by flow at right angles (in 
the front-to-back direction in the crucible). Perhaps more interesting, 
however, was the observation that the lateral flow produced by the 
lateral temperature gradient extended for a considerable distance into 
the flat temperature region. This would suggest that flat temperature 
profiles in a melt should extend well beyond the actual region of web 
growth to avoid significant flow under the growing tips of the bounding 
dendrites. 

2.3 Web Quality Analysis 

The primary goal of the material evaluation activity in the web 
growth program was to relate the web quality, primarily residual stress 
and dislocation density distribution, to the growth configuration and 
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growth conditions. Web structure evaluation techniques are compiled in 
Table 8. The technique for experimental evaluation of residual stress 
in a web crystal is depicted in Figure 19. 



Table 8. Web Structure 

Evaluation/Techniques 


PARAMETER 

MEASUREMENT TECHNIQUE 

1. 

Residual Stress 

Web Split Width Measurements 

2. 

Dislocation Density 

Etch Pit Counting 

3. 

Defect Type, Distribution 
and Structure 

X-Ray Topography 





The initial activity on this task under this contract involved 
organizing existing etch pit density (EPD) data from crystals grown in 
the two principal growth configurations being used for production and 
development: the J435 and the J460L. 

Residual stress measurement data for these configurations are 
presented in Figure 20. Figure 21 presents similar stress measurements 
along with a comparison of the average etch pit density as a function of 
position along the length of the web for two crystals, each from the 
J435 and the J460L configurations. Here it can be seen that in general 
the EPD for the J435 configuration is four to ten times larger than the 
EPD for the J460L configuration. Further, the EPD for the J435 material 
increases along the length of the crystal. In one of the J460L 
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Figure 20. Measured residual stress as a function of length along 

various crystals grown in J460L and J435 configurations. 


33 




J460 LID CONFIGURATION 



20 

Etch Pit 

IS 

Density 

12 


s 



20 

Residual 

10 

Stress 


CUdyn/cm*) 

0 



200 900 400 

Position, cm 


J435 LID CONFIGURATION 



40 

Etch Pit 

SO 

Density 
(10 3cm"*) 

20 


10 


0 


SO 

Residual 

Stress 

40 

(Udyn/cm*) 

20 


Z025-3 

2248-3 



Figure 21. 


O 100 200 900 400 SOO 000 

Position, cm 

Comparison of etch pit density and residual stress data from 
two web growth configurations. 







crystals, there is an initial peak in the EPD; however, this is followed 
by a decrease to a relatively constant level. 

The relationship between the EPD and the residual stress can be 
seen from data presented in Figure 21. The increase with length of the 
EPD for the J435 material is reflected in the residual stress data for 
the material, although the experimental scatter would make a strictly 
functional relationship rather vague. In the case of the J460L 
material, the initial peak of the EPD is echoed by an initial peak in 
the residual stress value; the relatively small overall value of the EPD 
is consistent with the small values of the residual stress. 

Topographic studies of specimens from this material show that 
EPD and stress values are generally consistent with the slip systems 
visible in the specimens. Additionally, there are the typical 
dislocation lines running along the length of the samples. The effect 
of these dislocations on stress and on the electrical properties of the 
material is still not known. 

In late September and October of 1984, crystal sampling from the 
development furnaces was begun on a routine basis. The initial sampling 
plan used is listed in Table 9. The EPD and stress data collected for 
the month of October are given in Table 10. It is important to note 
that the J and R furnaces had long susceptor systems with sloped wall 
crucibles, while the N and Z furnaces had normal length susceptors with 
sloped wall crucibles (all susceptors and work coils were rectangular). 
Also, all the runs were continuously replenished during crystal growth 
to maintain constant melt level. All the lids and shields were 
variations of the basic J460 design which were designed to produce low- 
stress material. 

With some exceptions, for example Run R470, the data in Table 10 
indicated that the average EPD in the most recent runs was larger than 
observed in crystals grown from earlier susceptor designs. In addition, 
the residual stresses were also larger and tended to be negative. Based 
on these limited initial data, the results tended to indicate that there 
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Table 9. Crystal Sampling Protocol 


Sample Selection: From all crystals grown greater than 3 m long: 

a. Button end of crystal 

b. 33 cm section from end of each succeeding meter 

c. Any other samples as determined by engineer 

Data Requirements: For samples provided as above, provide: 

a* Total Crystal Length 
b* Growth Velocity 

c. Melt Undercooling 

d. Melt Height 

e. Work Coll Settings 


had been a change in the temperature distribution in the lids and 
shields as the result of the use of rectangular susceptors and sloped 
wall crucibles. Note, however, that the stresses and EPDs were still 
much smaller than similar data for the J435 crystals. 

During the months of October and November of 1984, all the 

development furnaces were converted to sloped-wall systems of some 

nature; at the end of the period, only the N furnace had a short 

susceptor. The focus of the web quality analysis in November was to 

assess whether these system changes had any significant effects on the 

quality of the ribbon crystals. Initially, the etch pit densities and 

stresses for the J460L configuration in the long susceptor system were 

larger than for the same configuration using the short susceptor. 

Adjustment of the system temperatures, for example by coil position, 

brought the residual stress and etch pit densities back to values 

2 

comparable to the short susceptor cases, e.g., ±10 Mdyn/cm and less 
2 

than 6000 pits /cm as shown in Table 11. 

A systematic evaluation of etch pit density as a function of 
ribbon length was made for a number of long crystals grown in either the 
R or Z furnaces. Prior to adjustment of the system temperature, for 
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Table 10. Stress and Etch Pit Density Data From Web Grown 
with J460 Design Variations (October Data) 


Run-Crystal-Sample 

J535-1-1 

J535-1-4 

J535-1-6 

J535-1-7 

J539-1-3 

J539-2-3 

J539-8-B 

J539-8-4 

N113-2-A 

N1 13-2-5 

N117-1-B 

N117-1-4 

Z045-5-5 

R468-10-A 

R468-10-3 

R469-2-2 

R469-2-5 

R470-1-3 

R470-4-1 

R470-4-4 

R470-4-7 

R470-7-1 

R470-7-5 

R470-9-5 


Etch Pit Density 

1.2 X 10^ cm"^ 
5.8 

10.9 
114.0 

18.4 

23.9 
0.55 

12.5 
0.7 

3.7 

6.6 

12.5 

1.1 

0.3 

5.0 

5.6 

1.8 
0.7 
0.6 
0.9 

4.6 
1.5 

10.4 

3.0 


Stress 

2.3 Mdyn/cm^ 

9.5 

-36.9 

- 0.2 

-13.0 

-46.9 

-5.7 

-4.7 

6.0 

-9.7 

2.5 
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Table 11. Etch Pit Densities and Residual Stress from Web Grown 
After System Adjustments (November Data) 


Crystal . 
Run No. Piece Wo_. 


Residual Stress 
EPD xIO^ fMdvn/cm^) 


R461 

1.3 

R465 

1.8 

1.3 

1.7 

R407 

3.1 

3.9 

3.11 

R473 

1.3 

1.5 

R473 

6.2 

6.7 

6.11 

R474 

4.D 

4.3 

R474 

5.1 

R474 

6.1 

6.4 

R474 

8.9 

R474 

10. A 
10.1 
10.3 
10.7 

R475 

1.1 

1.6 

2047 

1.1 

1.6 

2048 

3.1 

3.8 

3.14 

2049 

4.1 

4.5 

4.8 

4.11 

4.14 

4.17 

2049 

7.16 

N117 

17.2 

N118 

1.1 

N119 

2.1 

2.5 

J539 

1.3 

J539 

2.1 


1.5 


.825 

4.33 

3.79 

-23.41 

.917 

18.6 

52.2 


2.76 

8.6 


3.2 

1.11 

1.5 


1.82 

.572 

8.54 

2.2. 


4.0 

31.9 

12.0 


14.7 

.628 

.458 

.484 

.471 

7.12 

.275 

12.9 

1.57 


4.03 

10.9 + 44.5 

4.99 

1.44 

2.54 


1.415 


.812 

3.09 

.903 


1.073 


1.336 


.864 

15.18 


7.79 

4.03 


7.7 

-8.48 


43.08 

5.62 


18.4 


4.52 
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2 

example in run R-467, etch pit densities exceeded 15,000 pits/cm and 
increased along the ribbon length. From later runs, for example runs 
R-470, R-473, R-474, Z-048, and Z-049, etch pit densities are small and 
essentially constant along the ribbon length (Figure 22). For several 

2 

of the ribbons, etch pit densities range from only 500 to 1500 pits/cm 
along many meters of crystals* This suggested that with the proper 
thermal configuration of the web growth system, long lengths of high- 
quality web crystals could be grown. 

Average stress and etch pit density data measured on crystals 
grown in sequential runs on the R furnace are presented in Figure 23. 
These runs were conducted immediately after introducing the long (10**) 
susceptor into the furnace. The data show a temporary increase in defect 
density initially after the modification, but the defect densities were 
substantially reduced by incorporating various system adjustments. 

2. A Experimental Web Growth 

The objectives of the experimental effort conducted with the 
developmental furnaces were to investigate high throughput growth using 
advanced furnace concepts and to demonstrate the ability to grow long 
crystals while operating in a steady-state replenishment mode. 

This task included the modification and development of existing 
web growth furnaces to attain higher area throughput rates. Also, 
improved melt replenishment techniques were devised to attain long-term 
continuous operation. Growth parameters were optimized by conducting 
tests and analyzing experimental growth characteristics. Figure 2A 
shows the experimental variations that are addressed for a typical 
growth configuration. 

The initial phase of the experimental growth program was devoted 
to converting all four developmental furnaces to a rectangular sloped- 
wall crucible configuration. Work performed under a previous contract 
indicated that better temperature stability and improved isolation of 
the growth and replenishment region of the melt are achieved using the 
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Figure 23. Average stress and etch pit density measured for R-furnace crystals. 
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Figure 24. Experimental variations addressed in typical growth hardware 
modifications. 



rectangular sloped-wall configuration. Work performed under this 
contract verified the improvement. 

In September 1984, 16 growth runs were made in the four primary 
web growth development furnaces, all with continuous melt replenishment 
A growth run consists of furnace set-up comprising many crystal starts 
and lasting up to a full week of operation. 

The main effort during September was aimed at obtaining 
acceptable profiles in the melt for the two sizes of rectangular growth 
systems. As might be expected, the susceptor shielding requirements 
were found to be somewhat different from those in the oval growth 
systems previously used for web growth. 

All of the growth configurations used during September were 
variations of the J460L. Principal variations included changes in the 
lid growth slot region and in the top shield spacing. These variations 
were introduced primarily to reduce the y and x stresses in the growing 
web* In addition, a vertical thermal element was added in one run in 
the N furnace, yielding some improvement in buckling stress over the 
baseline configuration, although the results have not yet been fully 
evaluated. 

Systematic run-to-run changes were also made in susceptor 
shielding to improve the lateral temperature distribution in the melt. 
Post-run examination showed that in all cases the sloped-wall crucible 
was in tight contact with the susceptor cavity, so that a significant 
improvement in the uniformity and reproducibility of heat transfer to 
the melt was achieved by conversion to the rectangular sloped-wall 
growth systems. 

The first three runs in the R furnace with the new extended 
growth system produced very encouraging results. A total of 70 meters 
of crystals over 1 meter long were grown with continuous melt replenish 
raent. Widths reached 4.8 cm. Growth problems attributable to melt 
replenishment were significantly reduced from previous experience. The 
purpose of these runs was to systematically move from a J460L 
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configuration to the lower stress J483L configuration, while making 
appropriate adjustments in the susceptor shielding to preserve an 
acceptable temperature distribution in the melt* The lid growth slot 
was sufficiently wide to allow the crystal to reach a width at which 
buckling would occur with the J460L configuration. Thus, the chief 
cause of crystal termination was the onset of stress induced by crystal 
degradation. Results showed an expected increase in the width at which 
deformation occurred. 

In all but one case, modifications in growth configurations gave 
the desired results or at least moved in the desired direction. One 
exception, however, was the J460L-M4 modification, which performed 
slightly poorer than the baseline J460L in terras of the buckling 
stresses. A further modification, designated the J460L-M5, was defined 
as a configuration lying roughly between the baseline J460L and the M4 
version. 


In October, 17 runs were made in the four development furnaces. 
Two types of lid configurations were run in the long rectangular 
configurations, the J460/J483LS designed for 5 cm crystals and the 
J460L-M3 designed for a 4 cm width limit. Variations in top shield 
configurations were tested for both lid configurations. Runs with the 
short rectangular systems were made with the J460L M3 and M5 
configurations with some variations in susceptor shielding to improve 
melt profiles. 

One run in the R furnace using the long J460/J483LS configura- 
tion with an expanded top shield spacing produced record results: the 

longest duration of continuous replenishment on a single crystal — 10.8 
hours; the longest replenished crystal — 6.8 meters; and, up to that 
time, the largest single-run replenished output — 36 meters. This run 
utilized a crucible with double barriers of a new design. 

Another run was made in the J furnace using the same lids but 
Incorporating a 2.5 cm vertical element in place of expanded top shield 
spacing to achieve the desired temperature gradients in the web. This 
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configuration also seemed to grow well, yielding a 4.7 meter replenished 
crystal. However, when attempts were made to combine both the expanded 
shield spacing and the vertical element, a configuration which should 
produce minimum stress in the web, the oxide deposition changed 
substantially, prohibiting growth. The convective gas flow patterns 
seem to be quite sensitive to the design of the vertical shield stack. 
Based on these results, work was initiated on designs to incorporate the 
desired web temperature effects, without imposing an unacceptable 
convective gas flow pattern. 

Growth in the short crucible systems was variable. A problem 
with front-back asymmetry was identified, and this condition was 
improved by changing the relative front-back heat losses with 
perforations in the susceptor side shields. Another sporadically 
encountered problem was welding of the adjustable end shields to the 
susceptor, so that they were no longer adjustable. Several methods for 
preventing shield contact with the susceptor were tested, and the most 
reliable appeared to be the use of ceramic separators at the critical 
points. 

One of the better runs in the N furnace with a J460L M3 
configuration yielded 36 meters of crystals, with a maximum replenished 
crystal length of 5.8 meters. 

During the course of this work, some distinct advantages were 
found with the long rectangular susceptor/crucible geometry. First, 
acceptable melt temperature profiles and good growth were found to be 
achievable with the higher stack/lower stress designs generated from the 
thermal models. Secondly, a significant improvement in decoupling 
between growth and replenishment effects was achieved. Thus, it was 
found easier to concentrate on those parameters that directly affect the 
quality of crystal growth. For these reasons, installation of a third 
elongated growth system was initiated in the Z furnace. 

In November, the Z furnace was converted to the long rectangular 
sloped-wall growth configuration. A new lid design was fabricated for 
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use in this system. The design, designated the J460LS-M6, is the same 
as the J460 in the vertical direction but incorporates a slot geometry 
aimed at a crystal width limited to 3.3 cm, well below the buckling 
stress limits for the J460 design. The rationale for this design is to 
grow long crystals at higher growth velocities. 

During the first run, crystal lengths of 4.6 and 5.4 meters were 
achieved. In the next run, a new record replenished crystal, 7.06 
meters in length, was grown. In addition, a 2.9 meter long crystal was 
grown at a growth velocity of 2.1 cm/min. 

The performance of this growth system was sufficiently encouraging 
that conversion of the N furnace to a long rectangular susceptor with lids 
and shields to duplicate the Z furnace configuration was initiated. 

The R furnace was run in November with the J460LS-M3 configura- 
tion, aimed at 4 cm crystal widths. The maximum crystal length was only 
5.4 meters. Some effort was directed toward improving growth velocity 
by testing growth at various accurately measured melt levels and coil 
elevations. A number of crystals were grown at 1.5 to 1.6 cm/min, below 
the predicted capabilities of the growth geometry but still an 
improvement over past averages. 

A new cover plate was installed on the J furnace. The new cover 
Incorporates a wider chimney, laser ports positioned to allow melt level 
sensing in the long susceptor, repositioned end-shield positioning rods 
to permit unencumbered measurement of dendrite edge thickness, provision 
for improved susceptor temperature measurements, and provision for 
covering view ports to reduce oxide accumulation. Modified electronics 
for the melt level control system were also installed. The hardware for 
the new laser mounting was received, and installation was initiated. 

Results of thermal modeling, discussed in Section 2.1 of this 

report, indicated that the J460 growth configuration is capable of 

growing web with sufficiently low buckling stresses to sustain growth at 

widths in excess of 5 cm and at speeds up to 2 cm/min. This combination 

2 

yields area throughput greater than 10 cm /min. These thermal models 
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address temperature distributions in the growing web, but melt 
temperature distributions must satisfy the requirements for growth 
before the potential of any furnace configuration can be met. 

Figure 25 is a plot of long crystals (longer than 4 meters) 
grown in the replenishment mode in a nine-month period during 1984. No 
crystals longer than four meters were grown in this mode prior to this 
time. The plot shows the increasing frequency of long crystals since 
the introduction of the long (10") rectangular sloped-wall system. 

Figure 26 is a pie-chart that identifies the types of termina- 
tions associated with the long crystals. The average growth time for 
these crystals ranging in length from 4 to 7 meters was 7.1 hours and 
the maximum was 11 hours. From this data it is obvious that melt 
replenishment problems have been significantly reduced through the 
introduction of the long rectangular, sloped-wall growth system. 

Figure 27 addresses "effective" growth velocities in develop- 
mental furnaces along with representative data from the pre-pilot line 
furnaces using the standard J435 configuration with no replenishment. 
Effective growth velocity is defined as the ratio of the total length of 
crystals grown in a given growth run to the total time that the furnace 
was held at operating temperature. A typical growth run lasts 
approximately one week. 

The data show that despite the experimental nature of the 
developmental growth runs, the effective growth velocity is nearly twice 
that of the standard pre-pilot line furnaces. The primary reasons for 
the low effective velocity in the pre-pilot line furnaces is the lost 
growth time involved with manually replenishing silicon between crystals 
and the necessity for gradually reducing the growth velocity as the 
unreplenished melt level falls. 

Table 12 summarizes the equipment modifications made to experi- 
mental growth furnaces during the course of this three-month program. 
Table 13 summarizes the lid and shield configuration variations investi- 
gated on the contract to improve stress and growth characteristics. 
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Figure 25. Record of crystals grown longer than 4 meters during replenishment. 
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Figure 26. Causes for termination of long crystals. 

Figure 28 shows the improvements made in terms of crystal length 
and replenished growth time since the initiation of work described in 
this report. Maximum crystal lengths have more than tripled, and growth 
times have increased by nearly a factor of four. Despite this progress, 
however, the 10 meter long crystal goal has not yet been reached. A 
primary factor in the lower than anticipated crystal lengths was the 
fact that growth speed improvements did not keep pace with growth 
stability improvements. It has been observed that the switch to the 
long crucible configuration that produced the improved growth stability 
(better temperature distribution in the molten silicon) had an adverse 
effect on axial temperature distribution along the length of the grown 
crystal; this produced higher stresses in the growing ribbons. To 
compensate for the increased stresses, it was found necessary to grow 
thicker web, hence lower pull speeds were used. 
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Effective growth velocities tor pre-pilot line and experimental growth fnrnaces 



Table 12. Major Equipment Modifications 


• Z furnace converted to long rectangular growth system in November 

e N furnace converted to long rectangular growth system in December 

• New top cover installed on J furnace with the following features: 

1 . Wider chimney 

2. Laser ports matched to long susceptor 

3. Allows for dendrite thickness monitoring 

4. Provisions for covering view ports 

5. Provisions for more temperature measurements 

o Laser melt sensing hardware Installed on J furnace 

• Feeder with commercial pellet material installed on R furnace 

• Larger susceptor concepts Identified 

1. Material requirements identified and procurement Initiated 

2. Equipment requirements identified and procurement initiated 


Table 13. J460L Type Lid and Shield Configurations 


Designation 


Design Function 


J460L 

J460L-M1 

J460L-M2 

J460L-M3 

J460L-M4 

J460L-M5 

J460LS 

J460LS-M3 

J460/483 LS 

J460LS M6 

NOTES: 1. 

2 . 


Low Stress Web at Controlled Width of ~4 cm 

Low Stress Web at Controlled Width of ~3.8 cm 

Alter Temperature Distribution in Y-Direction 

Alter Temperature Distribution in Y-Direction 

Rectangular Version of J460L-M2 

Modification of J460L-M2 

Stretched Version of J460L, 5 cm Wide Web 

Long System, M3 Lids 

Long System, J483 Top Shield Stack 

Long System, Controlled Width of 3.3 cm 

Mx designates lid slot modification 

S = Long Growth System = 10 Inch Crucible 
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Figure 28. Progress (crystal length and time) made toward 1984 
growth objectives. 
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3. RECOMMENDATIONS 


Work performed on this contract has led to a significantly 
Improved understanding of the parameters controlling the rate at which 
dendritic web silicon ribbon crystals can be grown. Realizing the 
improvements required to demonstrate commercial technical readiness of 
dendritic web growth requires additional effort in all areas covered by 
this contract. The substance of the recommended effort is described in 
this section of the report. 

3, 1 Stress Modeling 

Computer modeling must be conducted to define growth 
configurations capable of production of high-quality ribbon (having 
suitable thickness and low stress) at Increased area growth rate and 
under conditions of continuous melt replenishment. "Continuous melt 
replenishment" is defined as the condition in which the silicon melt is 
being replenished at the same rate that it is being consumed by ribbon 
growth, providing constant melt level (CML) conditions. 

The computer modeling should be directed at assessing web stress 
tendency of the web to buckle near the top of the shield stack (far 
stress) and plastic stresses generated very near to the growth interface 
(near stress). 

Modeling of stresses generated in growing web crystals near the 
crystal-liquid interface (near stress) should be performed to 
characterize and minimize plastic flow and plastic-elastic stress 
interactions. The coupling of in-plane plastic deformation, elastic and 
plastic buckling, and resulting web termination modes should also be 
investigated. This effort should include an assessment of available 
high-temperature mechanical property data, allowing expansion of the 
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existing WECAN model to include elastic anisotropy of silicon and yield 
stress criteria. 

3.2 Fluid Flow Modeling 

A study of the effects of convection in the molten silicon must 
be conducted to achieve the improved growth stability required for high 
speed growth. Both two- and three-dimensional analyses of the thermal 
environment and convective flow in the melt are required. Susceptor 
modifications to reduce temperature gradients producing convection 
should be assessed. For example, the effectiveness of incorporating 
baffles in the crucible to Inhibit melt convective currents should be 
evaluated. Finally, an analysis of the potential for magnetic 
stabilization of the molten silicon against convection should be made. 

3.3 Web Quality Analysis 

Planar and cross-section specimens of web using Transmission 
Electron Microscopy (TEM) should be examined to assess possible 
correlations of (111) surface etch features with the internal defect 
structure induced by deformation under web growth conditions. Growth- 
induced stresses must be related to both fine and large-scale defect 
structures which can degrade the web and minimize conditions that favor 
deformation. 

Features such as slip traces and location of contaminants (Si02» 
metals, etc.) that may result from gas phase interactions should be 
correlated with growth configurations. Techniques for minimizing those 
defects that adversely affect cell performance and influence plastic 
flow must be identified. 

3.4 Closed Loop Furnace Operation 

A closed loop control system must be devised and developed to 
allow long-term, stable, constant melt level growth of dendritic web 
crystals. This system must control the (1) average melt temperature. 
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(2) left-right melt thermal symmetry, and (3) melt level to provide 
steady state growth (no terminations due to pull-outs or extra 
dendrites). The system must maintain constant melt level either with an 
active laser-sensor control and silicon pellet feeder or by fixing the 
feed rate of a calibrated pellet feeder to exactly match the silicon 
consumption during crystal growth. In conjunction with this effort, 
alternate methods for pellet feeding, sensing of pellet drops (uniform 
dispensing), and improved crucibles and lid configurations to maintain 
thermal symmetry while feeding should be pursued. Alternate pellet- 
feeding methods investigated should include vibratory techniques. 
Acoustic and photoelectric sensors should be investigated to detect and 
count pellet drops. 
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